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Q Abstract 

The oversampling multiscale finite element method (MsFEM) is one of the most popular methods 
for simulating composite materials and flows in porous media which may have many scales. But the 
method may be inapplicable or inefficient in some portions of the computational domain, e.g., near 
the domain boundary or near long narrow channels inside the domain due to the lack of permeability 

i— i information outside of the domain or the fact that the high-conductivity features cannot be localized 

within a coarse-grid block. In this paper we develop a combined finite element and multiscale finite 

y-^- element method (FE-MsFEM) , which deals with such portions by using the standard finite element 

method on a fine mesh and the other portions by the oversampling MsFEM. The transmission con- 
ditions across the FE-MSFE interface is treated by the penalty technique. A rigorous convergence 
analysis for this special FE-MsFEM is given under the assumption that the diffusion coefficient is 

aperiodic. Numerical experiments are carried out for the elliptic equations with periodic and ran- 
i dom highly oscillating coefficients, as well as multiscale problems with high contrast channels, to 

demonstrate the accuracy and efficiency of the proposed method. 
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1 Introduction 

Let f2cM ri ,n = 2,3bea polyhedral domain, and consider the following elliptic equation 

-V • (a e (x)Vu e (i)) = f(x) in Q, 



u e (x) = ondil, 



where < e 1 is a parameter that represents the ratio of the smallest and largest scales 
in the problem, and a t (x) = (a^-(x)) is a symmetric, positive definite, bounded tensor: 

(1.2) \\e < a'p-X.C, < AKI 2 V£ E R n , x E ti 

for some positive constants A and A. 
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Problems of the type (1.1 ) are often used to describe the models arising from composite 
materials and flows in porous media, which contain many spatial scales. Solving these 
problems numerically is difficult because of that resolving the smallest scale in problems 
usually requires very fine meshes and hence tremendous amount of computer memory 
and CPU time. To overcome this difficulty, many methods have been designed to solve 
the problem on meshes that are coarser than the scale of oscillations. One of the most 
popular methods is the multiscale finite element method (MsFEM) [2E1 1121 S3], which 
takes its origin from the work of Babuska and Osborn [HI [8]. Two main ingredients of the 
MsFEM are the global formulation of the method such as various finite element methods 
and the construction of basis functions. The special basis functions which constructed 
from the local solutions of the elliptic operator contain the small scale information within 
each element. By solving the problem (1.1) in the special basis function space, they get a 
good approximation of the full fine scale solution. We remark that there are many other 
methods proposed to solve this type of multiscale problems in the past several decades. 
See, for instance, wavelet homogenization techniques [HI [30], multigrid numerical ho- 
mogenization techniques (35] [50], the subgrid upscaling method [21 E], the heterogeneous 
multiscale method [211 E21 [23] , the residual-free bubble method (or the variational multi- 
scale method, discontinuous enrichment method) [131 ED H23 Ell 051 [51] , mortar multiscale 
methods [U [53] , and upscaling or numerical homogenization method [201 E21 [59] . We re- 
fer the reader to the book [26j for an overview and more other references of multiscale 
numerical methods in the literature, especially a description of some intrinsic connections 
between most of these methods. 

In this paper, we focus on the MsFEM. Many developments and extensions of the 
MsFEM have been done in the past ten years. See for example, the mixed MsFEM 
dUE], the MsFEMs for nonlinear problems [271 [29], the Petro-Galerkin MsFEM [B], 
the MsFEMs using limited global information [251 EI]) an d the multiscale finite volume 
method [46J. In [43J, it is shown that there is a resonance error between the grid scale and 
the scales of the continuous problem. Especially, for the two-scale problem, the resonance 
error manifests as a ratio between the wavelength of the small scale oscillation and the 
grid size; the error becomes large when the two scales are close. The scale resonance is 
a fundamental difficulty caused by the mismatch between the local construction of the 
multiscale basis functions and the global nature of the elliptic problems. This mismatch 
between the local solution and the global solution produces a thin boundary layer in 
the first order corrector of the local solution. To overcome the difficulty due to the 
scale resonance, an oversampling technique was proposed in [121 [2H]. The basic idea is 
computing the local problem in the domain with size larger than the mesh size H and 
use only the interior sampled information to construct the basis functions. By doing this, 
the influence of the boundary layer in the larger domain on the basis functions is greatly 
reduced. 

However, for the coarse-gird elements near the boundary, in order to construct the 
multiscale basis functions, the oversamping MsFEM needs to assume that there is enough 
information available outside of the research domain, which is not applicable in practice. 
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To handle this problem, the natural way is to use the standard multiscale basis functions 
instead of the oversampling multiscale basis functions in the coarse-grid elements adjacent 
to the boundary, hence in this area we don't need to use the information outside the 
domain. We call this method as the mixed basis MsFEM. Since in the elements near the 
boundary, we use the multiscale basis functions without oversampling technique, the scale 
resonance comes out again hence pollute the accuracy. 

To overcome this difficulty, we introduce a new method in this paper which can improve 
the accuracy significantly. The proposed method separates the research area into two sub- 
domains such that one of them is contained inside the domain with a distance away from 
the boundary. Then in the interior sub-domain the oversampling multiscale basis functions 
on a coarse mesh (with mesh size H) are used. While, in the other sub-domain which 
adjacent to the boundary the traditional linear FEM basis functions are used on a mesh 
(with mesh size h) which is fine enough to resolving multiscale features. The difficulty to 
realize this idea is how to joint the two methods together without losing accuracies of both 
methods, i.e., how to deal with the transmission condition on the interface between coarse 
and fine meshes efficiently Thanks to the penalty techniques used in the interior penalty 
discontinuous (or continuous) Galerkin methods originated in 1970s [T0" | fTT | [l~9l 1561 l5| 16]. 
we may deal with the transmission condition on the interface by penalizing the jumps 
from the function values as well as the fluxes of the finite element solution on the fine 
mesh to those of the oversampling multiscale finite element solution on the coarse mesh. 
A rigorous and careful analysis is given for the elliptic equation with periodic diffusion 
coefficient to show that the if terror of our new method is just the sum of interpolation 
errors of both methods plus an error term of O(^) introduced by the penalty terms, 
where H is the mesh size of the coarse mesh. We would like to remark that besides the 
applications of penalty technique to the interior penalty discontinuous (or continuous) 
Galerkin methods, this technique is also applied to the Helmholtz equation with high 
wave number to reduce the pollution error [331 EH EZ1 ED] and applied to the interface 
problems to construct high order unfitted mesh methods [48, 58J. 

The other potential application of our proposed method is to solve the multiscale prob- 
lems which may have some singularities. For example, the multiscale problem with Dirac 
function singularities, which stems from the simulation of steady flow transport through 
highly heterogeneous porous media driven by extraction wells [16J, or the multiscale prob- 
lems with high- conductivity channels that connect the boundaries of coarse-grid blocks 
[38| l39~t I2H |52]. Our new FE-MsFEM may solve such problems by using the traditional 
FEM on a fine mesh near the singularities (and, of course, near the domain boundary) 
and using the oversampling MsFEM in the other part of the domain. To demonstrate 
the performance of the FE-MsFEM, we try to simulate multiscale elliptic problems which 
have fine and long-ranged high-conductivity channels. We remark that this kind of high- 
conductivity features cannot be localized within a coarse-grid block, hence it is difficult to 
be handled with standard or oversampling multiscale basis. The numerical results show 
that the introduced FE-MsFEM can solve the high contrast multiscale elliptic problems 
efficiently. The convergence analysis for multiscale problems with singularities and ap- 



3 



plications of the proposed FE-MsFEM to practical problems such as two-phase flows in 
porous media and other types of equations are currently under study. 

The rest of this paper is organized as follows. In Section |2j we formulate the FE- 
MsFEM for the model problem. In Section [3j we review some classical homogenization 
results for the elliptic problems and give an interior H 2 norm error estimate between the 
multiscale solution and the homogenized solution with first order corrector. In Section [4], 
we give some approximation properties for the oversampling MsFE space and the linear 
FE space, respectively. The H l error estimate of the introduced FE-MsFEM is given 
in Section [5j In Section [6] we first give some numerical examples for both periodic and 
randomly generated coefficients to demonstrate the accuracy the proposed method, and 
then apply our method to multiscale elliptic problems which have fine and long-ranged 
high-conductivity channels to demonstrate the efficiency of the method. Conclusions are 
drawn in the last section. 

Before leaving this section, we fix some notations and conventions to be used in this 
paper. In the following, the Einstein summation convention is used: summation is taken 
over repeated indices. L 2 (Q) denotes the space of square integrable functions defined in 
domain Q. We use the L 2 (Q) based Sobolev spaces H k (Q) equipped with norms and 
seminorms given by: 




IM|wk.«>(fi) (l M lw fc '°°(n)) is the W k '°° norm (seminorm) of umfl. Throughout, C,C\,C2, ■ ■ ■ 
denote generic constants, which are independent of e, H and h unless otherwise stated. 
We also use the shorthand notation A < B and B > A for the inequality A < CB and 
B > CA. The notation A ~ B is equivalent to the statement A < B and B < A. 

2 FE-MsFEM Formulation 

In this section we present our FE-MsFEM. We describe the method only for the case of 
dealing with the difficulty of lack information outside the domain in the oversampling 
MsFEM. Of course, the formulation can easily be extended to the case of dealing with 
singularities. 

We first separate the research area Q into two sub-domains Qi and Q2 such that 
02 CC n and Q = Q\ U Q2 U T, where T = dQi fl 8^2 is the interface of Qi and Q2 
(cf. Fig. [T]). For simplicity, we assume that the length/area of T satisfies |T| = 0(1). 
Let M. h be a triangulation of the domain Qi and M. H be a triangulation of the domain 
Vt 2 , and denote and V ^ the two partitions of the interface V induced by Aih and 
M.H) respectively. We assume that on the interface T, A4 H and M.^ satisfy the matching 
condition that Th is a refinement of Th- Clearly, each edge/face in Th is composed of 
some edges/faces in Th- Combining the two triangulations together, we define Aih,H as 
the triangulation of f2 (See Fig. [T] for an illustration of triangulation). For any element 
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K G M.h (or K G M-h), w e define /i^ (or Hk) as diam(K). Similarly, for each edge/face 
e of K e G .M^ (or i£ of 1^ G ), define h e as diam(e) (or H E as diam(.E')). Denote 
by h = m&XKeMh h-K and H = max KeMlI H K . We assume that h < e < H, that {Alff} 
and {M.h\ are shape-regular and quasi-uniform. 




Figure 1: A separation of the domain and a sample mesh. 

For any point on T, we associate a unit normal n, which is oriented from Q\ to f^. We 
also define the jump [v] and average {v } of v on the interface T as 



(2.1) [v] := v|m - u|n 2 , 

Introduce the "energy" space 



(2.2) 1/ := {f : v\n t = v u where t>; G #o(fi) H # s (fi), i = 1,2}, for some s > 



Testing the elliptic problem ( 1.1[ ) by any v £ V, using integration by parts, and using the 
identity [vw] = {v} [w] + [v] {w}, we obtain 

2 

[ a e VM e ■ Vv - / {a f V« e • n} [v] = / fv. 

i=1 JUi Jr Jn 



■5 



Define the bilinear form Ar{-, ■) on V x V: 



(2.3) A p (u,v):= / a e Vw-Vi; 

KeM h$H K 

(2.4) ~J2 J { ' n ^ M + ^ N ' n ^ ) 



+ J (u,v) + Ji(u,v), 

(2.5) J (u,v):=J2^ [ MM, 

(2.6) Ji(m, u) := ^ 7l p / [a £ Vw • n] [afVv ■ n] , 

eer\ ^ 

where /3 is a real number such as —1, 0, 1, and 70, 71, p > will be specified later. Define 
further the linear form F(-) on V: 

F{v) := / fv. 



It is easy to check that the solution u e to the problem (1.1 ) satisfies the following formu- 
lation: 

(2.7) A p (u e ,v) = F(v) VveV. 

To formulate the FE-MsFEM, we need the oversampling MsFE space on A4h defined 
as follows (cf. |15l HH 125]). For an y K e M H with nodes {xf}^, let {(ff}^ be 
the basis of P\{K) satisfying ipf(xf) = 1 < i,j < n + 1, where 5y stands for the 
Kronecker's symbol. For any K e A^h, we denote by 5 = S^iT) a macro-element 
(simplex) which contains K and satisfies that H s < C\H K and dist(K, dS) > C H K , 
where C% > is independent of Hk and Co is from (2.17). We assume that the macro- 
elements S(K) are also shape- regular. Denote by {<pf the nodal basis of Pi(S) such 
that ipf(xj) = Sij, 1 < z, j < n + 1, where are vertices of S. 

Let e H l (S),i = 1, • • • , n + 1, be the solution of the problem 

(2.8) _V-(a e V^f)=0 in S, ^ \ as = (pf . 

The oversampling multiscale finite element basis functions over .fT is defined by 

(2.9) fa* = c$ipf\ K in K, 
with the constants so chosen that 

(2.10) V f = c%^\ K in K. 

The existence of the constants cfj is guaranteed because {(pj}™^ 1 also forms a basis of 
P 1 (K). 
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Let OMS(-ff) = span {ipi be the set of space functions on K. Define the projection 
U K : OMS(il~) -> P X {K) as 

u K ip = dtp? if v = e OMS(K). 

Introduce the space of discontinuous piecewise "OMS" functions and the space of discon- 
tinuous piecewise linear functions: 

X H = \jpH : Mk e OMS(K) ViT G M H }, 
W H = {w H : w h \k G P 1 (K) \/K G M H }. 

Define Uh '■ Xh — > Wh through the relation 

(2.11) U H ip H \ K = U K ip H for any K G M.h, £ Xh- 

The oversampling multiscale finite element space on is then defined as 

X H = {ipH e X H : U H ip H G Wh C i? 1 ^)}, 

where W# = D ii^^) is the if ^conforming linear finite element space over A4h- In 
general, X H <£. H 1 (fl 2 ) and the requirement Hh^h G W^h is to impose certain continuity of 
the functions ipn £ -Xi? across the inter-element boundaries. According to the definition 
of 11^, we have IIk-^ = <£>f . Since <pf is continuous across the element, this above 
requirement Hh^h £ We is satisfied naturally since in each node we only have one 
freedom ( unknowns ). 

Denote by Wh the if ^conforming linear finite element space over M. h and by 

(2.12) W% := {w h EW h : w h = on dtt^T} . 
We define the FE-MsFE approximation space V^e as 

(2.13) V KH ■= {vh,H ■ v^hW = v h , v h}H \n 2 = v H , where v h G W%, v H G X H ] - 

Note that Vh,H V- We are now ready to define the FE-MsFEM inspired by the formu- 
lation (2.7): Find Uh,n G Vh,H such that 

(2.14) A fi (u h>H ,v hjH ) = F(v htH ) Vv h>H G V hjH - 

Remark 2.1. (a) If (5 = 1, then the the bilinear form Ap is symmetric and, as a conse- 
quence, the stiffness matrix is symmetric as well. If (5 ^ 1, e.g., (3 = —1, then the method 
is nonsymmetric. 

(b) The parameter p > satisfies that p < e. In fact, it is chosen as e in our later 
error analysis, while in practical computation, it may be chosen as the mesh size h. 

For further error analysis, we introduce several concepts related to the interface T and 
some discrete norms. Define the set of elements accompanying with the interface partition 
T h (or T H ) as following: 

(2.15) Kr h := {K G M h ■ K has an edge/face in T h }, 

(2.16) K Vh := {K G Mh ■ K has an edge/face in T H }. 
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Clearly, the number of elements in Kr h is 0(^=r) and the number of elements in Ky h is 
0( H n-i )- Moreover, we assume that 

(2.17) dist{r, dVt] > C H >h + 2e>0, 

where Co is a constant. Thus, we can define a narrow subdomain Qr CC Q surrounding 
r as 

(2.18) U r :=ru{i: x e Q±, dist(x, T) < h + 2e} 

U {x : x G fla, dist(x, r) < if + 2e}. 

Denote by 

(2.19) fir H = U{X : G ^r H }, fir h = U{-K" : K G tf r J. 

It is clear that Q-p H ,^r h C fir and dist(fir H ? <9fir) ; dist(Or fc , 9fir) > 2e, repectively. 
Denote by 

Nii lff =f e ii( ae ) i/2 ^nV)) v ^ n 

\K£Mh / KGA4 H 

We introduce the following energy norm and the broken norm on the space Vh,H- 

.2 \V2 
L 2 (fii) 



(2.20) (11^11^+ || (a e ) 1/2 V 

(2.21) lll^lll l h H := (\\v\\l + E ^ II M & (e ) + E 7i P II [a e Vv ■ n] \\% 

+ E^HKV,.n}||J 2(e) ) 1/2 . 



3 The homogenization results 

In this section, we assume that a e (x) has the form a(x/e). Moreover, we assume that 
aij(y) G Cp(M n ), where Cp(M. n ) stands for the collection of all C 1 (M n ) periodic functions 
with respect to the unit cube Y. It is shown that under these assumptions (cf. [121 E!)> 
u 6 converges in a suitable topology to the solution of the homogenized equation 



(3.1) 
where 



-V • (a*V«o(a:)) = f(x) infl, 
u (x) = ondQ, 



If ( dx j 

(3-2) a*j = — a ik {y) \ 5 kj + q^(v) ) d V- 
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Here is the periodic solution of the cell problem 

(3-3) -Vy(a(y)V yX J (y)) = Vy(a(y)e 3 ), j = l,---,n 

with zero mean, i.e., j Y X*dy = 0, and is the unit vector in the jth direction. The 
variational form of the problem (3.1) is to find Uo(x) G Hq(Q) such that 

(3.4) (a*ViZo,Vv) = (/,«) Vue^fi). 



It can be shown that a* is positive definite. Thus by Lax-Milgram lemma, (3.4) has a 
unique solution. If f(x) G L 2 (Q), from the regularity theory of elliptic equations, we have 



(3-5) \\ u o\\ H 2(n) ^ C \\J \\L 2 (n) ■ 

Let 9 e denote the boundary corrector which is the solution of 

-V • (a e V0 £ ) = in Q, 

(3-6) , ,duo(x) 

e = -X J {xe)— ondQ. 

dXj 

From the Maximum Principle, we have 
(3-7) lre|li°°(n) ~ l u olvi/i,cx>(Q) ■ 

In the following part, for convenience's sake, we will set 

du (x) 



(3.8) Ui(x, x/e) = Uq(x) + ex 3 (x/e)- 



dxj 



The following error estimates are known (cf. e.g. [14, 49, 16J). 

Theorem 3.1. Assume that uq G H 2 (Q) fl W 1,(X) (Q). Then there exists a constant C 
independent of e, the domain fl, and the function f such that 

(3.9) || V(u £ -u x - e9 e ) \\ L 2 {n) < Ce \u \ H2{n) . 
Moreover, the boundary corrector 8 e satisfies the estimate 

(3.10) ||eV0 6 ||i3(n) < Ce \u \ H 2 {n) + C^Je \3Vt\ |uol w i.~(n) ■ 
where \dVL\ stands for the measure of the boundary dfl. 

The following regularity estimate is an analogy of the classical interior estimate for 
elliptic equations in |40j Theorem 8.8, P. 183]. 

Lemma 3.2. Let w G H l (D) be the weak solution of the equation 

-V • (a e (x)Vw) = g(x), x G D, 
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where a e (x) satisfies (1.2) and g G L 2 (D). Then for any subdomain D' CC D, we have 
w G H 2 (D') and 

( 3 - n ) \ w \h*(d>) ^ [j, + I J \ w \m(D) + \\g\\ L i{D) 

where d' = dist(D', dD) and r is a constant such that |a e (x) — a e (?/)| < r \x — y\. 

Remark 3.1. The H 2 norm interior estimate for elliptic equations is well-known in the 
literature. The importance in above estimate is the explicit dependence of the bound on d' 
and r which is crucial in our analysis. 

Proof. First, we define the difference quotient as follow: 

A h v(x) = Aiv= V -^- +hek) - v{x) 



h 

where e&, 1 < k < n is the unit coordinate vector in the Xk direction. And then, following 
the proof presented in Theorem 8.8, P. 183], we obtain 



A 



/ |£V(A\i;)| 2 dx 

< (j. \NM\m D ) + hWmv)) (l|£V(A^)|| i2(D) + 2||AW£| 



L2(D) 



+ C\\iV(A h w)\\ L2{D) \\A h wV^\\ L2iD) 



where £ G C%(D) is the cut-off function such that < £ < 1, £ = 1 in D', and |Vf | < 2/df 
in D. By use of the Young's inequality and the Lemma 7.23 in |40j P. 168], it follows that 

\\Z Ah (Vw)\\ LHD) <c{^- \\Vw\\ L2{D) + |M| L2(D) + W^wViW^ 

- C (l + J') ^ W ^HD) + C \\9\\^(D) ■ 

By the L emma 7.24 in gO, P. 169] we obtain Vw G H 1 (D 1 ), so that w G H 2 (D') and the 
estimate (3.11) holds. This completes the proof. □ 

Utilizing the above H 2 interior estimate to equation — V-(a*ViZoay) = fxp j — lj ■ ■ ■ > n j 
where Xj is the coordinate variable in the jth direction, we obtain 



Lemma 3.3. Let u be the solution to (3.4). Assume that Uq G H 2 (D) and f G 
Then for any subdomain D' CC D with d' = dist(D', dD), we have 

(3-12) \ u o\h3(d>) ~ ^7 \ u o\h 2 (d) + II v /IIl 2 (d) • 



Further, by use of the H interior estimate (3.11 ), we obtain an H semi-norm interior 
estimate of the error u e — u\ in the narrow domain Q-p- 
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Theorem 3.4. Assume that up G H 2 (tt) n W 1,00 (n), and f\ Qr G H 1 ^). Then for any 
subdomain Q' CC fir with dist(Q', <9fir) > 2e ; we have 



(3.13) 



M il//2(n') ~ \ u o\mm + e II Wllz^r) + 



(«r) 



Proof. It is shown that, for any </j G i?o(fi) ( see [El p. 550] or [T5J p. 125]), 
(a(x/e)V(u e -wi),V^ n 



(3.14) 



a*V« , Vy?)n - a(x/e) V w + ex 



• du 
dx k 



aij(x/e)x k 



d 2 u dip 



— e 



a*Jx/e) 



d 2 u dip 
dxjdxk dxi 



dxjdxk dxi 

where oc^Ay) G Hl oc (R n ) are F-periodic and dependent only on the coefficients a.(y) (see 
[HI p.6]). According to the assumption, there exists a subdomain fi c CC O r such that 
fi' CC fi c with d = dist(fi', <9fi c ) > e and di = dist(fi c , <9fir) > e. From (3.14), it follows 
that in fi c , 



V • (a(x/e)V(u 6 - wi)) = e 



d_ 

dxi 



d 2 u 
dxjdxy 



otAx/e) 



d 2 u 
dxjdxk 



Thus, from Lemma 3.2, it follows that 



u e - Ui m 



< 



1 



dp 



\ue - ^i|jfi(n c ) + \ u o\h 2 (q c ) + e l M o|jf3(n o) 



Hence, from Theorem 3.1 and Lemma 3.3, it follows the result (3.13) immediately. □ 

We conclude this section with a local H 2 semi-norm estimate for u\ in fir h , which will 
be used in the convergence analysis. 

Lemma 3.5. Assume that Up G H 2 (Vt) n W X '°°(Q), and f\ Qr G H 1 ^). Then, 



(3.15) 



\ u i\m{si Vh ) £ KI H 2(n r ) + e 



Wo 



+ e||V/|| i2 



(n r ), 



where Qr h an d fir ore defined in (2.19) and (2.18), respectively. 
Proof. It is easy to see that 

;du 



\Ul\ H 2 iQrh) < \Up\ H 2 {nrh) 



+ 



ex 



< 



+ € 



dXj 
Our 



dxj 



+ 6 [wolira 



Thus, from Lemma [3 .3[ it follows that 

Mjjvnr, ) ~ l M o| H 2 (f7r) + e||V/|| L 2 (nr) + e- 1 \tt Fh \ 1/2 \\Vu \ 



L°°(n r .) ' 



which, combining with the fact that |firj = yields the result (3.15). 



□ 
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4 Approximation properties of the FE-MsFE space V^h 

In this section, we give some approximation properties for the oversampling MsFE space 
Xh and the linear FE space Wh, respectively. 

4.1 Approximation properties of oversampling MsFE space Xjj 

We first recall a stability estimate for IT# (see pj)J Lemma 9.8] or |28j Appendix B]). 

Lemma 4.1. There exist constants 7 and C independent of H and e such that if Hk < 7 
and e/H K < 7 for all K £ A4 H , the following estimates are valid 

CT 1 \\Vn. H w H \\ L 2^ K ) < \\Vw H \\ L 2 {K) < C \\VI1 h wh\\ L 2 {K ) Vvj h G X H . 

Next, we give some approximation properties of the oversampling MsFE space OMS(A). 

Lemma 4.2. For any K G M.h, there exists <pH G OMS(K) such that the following 
estimates hold: 

(4.1) \m - 4>h\ h i {K ) < H K \u \ H 2 {K) + eH^ 2 1 \uq\wi>°°(k) > 

(4.2) \\m - 4>h\\ L 2 {k) < H 2 K \u \ H2{K) + eH]/ 2 |«olwi.°°(K-) ' 

(4.3) \ui - 4>h\ H 2( K ) < e 1 H K \u Q \ H 2( K ) + H 7 ^ 2 1 Iwolwi,*^ + e \uq\ h z^ . 
Proof. We take 

(4.4) ct> H = Y, Mx?)Ux), 

node of K 

then 

Rk4>h = IhUo, 

where In '■ C(f2 2 ) — > Wh is the standard Lagrange interpolation operator over linear finite 
element space. By the asymptotic expansion, we know that 

(4.5) ( f>H = lHU + e X > d ^ + e6^ 



where 9~ G H 1 (S) is the boundary corrector given by 

(4.6) -V-(a e V0f) = O in 5, 0f| fl5 = - X j ^ 

By the Maximum Principle we have 

(4-7) 1^0 1 W X >°°(K) ' 
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which together with the interior estimate in Avellaneda and Lin [7J Lemma 16] imply that 



(4.8) 

Therefore 
(4.9) 



IL°°(K) ~ ||^e ||l°°(S) ~ xx K I " , U\W 1 >°°(K) ' 

jd(I H u ) 



V - Ihu - ex 



L\K) 



< H K |m q | w i, 



A' 



1/2 



< o-n/2-1 I I 



Further, since 



|V(W -^0)|| L 2 m <#A-| U o|jf2 



7 (9(m - J g M ) 



2 W 



< (H K + e) \u \ H 2 {K) 



we obtain the result (4.1) immediately. 



To prove the estimate (4.2), we first notice that, from (|4.7|), 



< Tjn/2 I I 



Therefore 
(4.10) 

Further, since 



\L*{K) 



h ~ IhUq + ex 



7 g(/ g M ) \ 

dxj J 



Uo\ 



L\K) 



ex 



\\u - I H u \\ L2{K) < H 2 K \u Q 
3 d(u - I H u 



L 2 (K) 



H 2 (K) ' 

< tH K \uq\ H 2^ , 



we obtain the result (4.2) immediately. 

To prove the estimate (4.3), it is easy to see that 

(4.11) 



(4.12) 



Uq - IhU \ H 2( K ) < \ u o\ H 2(K) > 

d(u - I h uq) 



ex 



dXj 



H 2 (K) 



< (1 + e 1 H K ) \u \ H2{K) + e \u Q \ HZ[K) . 



Let D be a simplex such that 

K CC D CC S and dist(A, dD), dist(D,dS) > H K . 



From Lemma 3.2 and following the proof of (4.8), we have 

\ H 2 {K) ^ {H K l + e _1 ) \dc\ H i {D) ^ h k 2 1 ( h k 1 + e_1 ) II ^ u o\\l°°(k) > 



(9 



which, combining with (4.11) and (4.12), yields the result immediately. This completes 
the proof. □ 
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From Lemma 4.2 we have the following local approximation estimates in Kr H - 
mn 

(4.13) 



Lemma 4.3. There exists ipn £ Xh such that, 

1/2 



(4.14) 
(4.15) 



( Yl \\ u i-M\v(iq) ^ H2 \ u v\m{n VH ) + e ^H\ u o\w^{n VH ) 
( 2 Y /2 e 

( k-V^I 2 ^ 



1/2 



< — |«o|H2(n r ) + 6 II V /llL2(Q r ) + \ U o\wi'™(n rH ) 



where Qr H is defined in (2.19). 
Proof. We take 

(4.16) i> H = 



pj node of Mh 



then ipH\ = 4>Bi which is defined in Lemma 



IK 



4.2 



From (4.2), it follows that 



22 IK - iI>h\\vhk) ^ H 4 \ u o\ 2 H2(n r „) + e2Rn l M o|^i,oo (n } Y X 



KeK Tf 



< £7 4 I 1 2 | 2 rr I 1 2 



which yields (4.13) immediately. Note here we have used the fact that the number of 

1 1 ipi . . . . 

elements in Kr H is Q{ ^n-i )- Similarly, (4.14) follows from (4.1). 

It remains to prove (|4.15). From (4.3) and Lemma 3.3, it follows that 



\ Ul ~ ^h\ 2 H 2 (k) < e 2 H 2 \u \ 2 H 2 [n ) + e 2 \u \ Hl 



+ 2 KlJl/l.oo^) 1 



K(LKy , 



< {e~ 2 H 2 + 1) \u \ 2 + e 2 \\Vf\\ 2 + H- 1 \u f wl , 



3 («r w ) 



□ 



which yields (4.15) immediately. The proof is completed. 

From Lemma 4J2 Theorem |3.1 , via taking the same ipn * n Lemma 4.3[ we also have the 
following result which gives an approximation estimate of the space X H (cf. [2H1 l2"6"t ITS]). 

Lemma 4.4. There exists ipn £ Xh such that, 



(4.17) \\u e - ip H \\ liH ^ H \u \ H 2 {n) + \u \ w i,oo { n 2) + \u \ w i : . 



■(0) 



14 



4.2 Approximation properties of linear FE space Wh 

Since |f2i|, the area/volume of Q,\, may be small, we prefer estimates with explicit de- 
pendence on it. To attain this aim, we use the Scott-Zhang interpolation instead of the 
standard FE interpolation in this subsection. 

We first introduce the Scott-Zhang interpolation operator Zh : H^ifli) — > W®, where 
H^(Qx) := {v G H 1 ^) : u = 0on dQi \ T}. For any node z in M h , let <f> z (x) be the 
nodal basis function associated with z and let e z be an edge/face with one vertex at z, 
then the Scott-Zhang interpolation operator is defined as |55j : 

(4.18) Z h v= ( I VdG^), 

nodezinMh ^ e ' A ' 

where ip z (x) is a linear function that satisfies J e ip z (x)w(x) = w(z) for any linear function 
w(x) on e z . Suppose e z C dfli for z G d£l\. It is easy to check that H^zll^^) < for , 

IIW*|L»(n,) ~ ^e"' Wllere : = SU PP(0 Z ) ; and 

(4.19) Z h v = v Vf G W£. 

This operator enjoys the following stability and interpolation estimates (see [55J): 
Lemma 4.5. For any K G Aih, we have 

( 4 -20) \\Z h v\\ Lx{K) < \\v\\ Lac(R) , \\VZ h v\\ Loc(K) < \\Vv\\ Loa{R) , 

(4-21) ||f - Z h v\\ L 2 [K) + h K \\v - Z h v\\ H1{K) < h\ \v\ H2[k) , 

(4-22) ||u - Z h v\\ Loo{K) < h\ \v\ w%00{R) , 

where K is the union of all elements in Aih having nonempty intersection with K. 

Moreover, we need the following error estimate between u\ and its Scott-Zhang inter- 
polant which uses only the regularity of the homogenization solution uq. 

Lemma 4.6. For any K G Aih, we have 



hU, , h K /2+1 



( 4 - 23 ) K - z h u i\ H i (K) < (h K + e + — ) \uq\ h2{ k) + 



\u 



where K is defined in Lemma 4-5 

dx 

ui - Z h ui\ m{K) < \u - Z h u \ H1{K) + e \x ] Vj - Z h (tfvj) \ H1 



Proof. Denote by Vj := f^. It is easy to see that 



(K) 



< 



\uq - Z h u \ H1{K) + e |(x J - Z h x 3 )vj\ H1{K) 

+ e \Z hX J (vj - (vj) K ) \ H1{K) + e \Z hX j (v J ) K - Z h (x j Vj) \ R1(K) 



T + II + III + IV. 
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where (•)# = rjU J k (') dx. From (4.21), we have 
(4.24) I < h K \u \ H 2 



(JO • 



From the facts that \Vx j \ < e -1 , |VV'I < e" 2 and ( |4.20[ )-( [422| ), we have 

n/2+1 



(4.25) 
(4.26) 



II <> 



h 2 



l u 0|iyi>°°(X) + ~ l M o|//2(^) 



III < (/ix + e) \u \ H 2 



(K) 



where we have used the Poincare inequality to derive the second inequality. It remains to 
estimate IV. According to the definition of Scott-Zhang interpolation, we have 



IV = e 



node zdK 

5, eh J- 



m(K) 



^( Z ) / fa* 1 i( V j)K- V j) 

W*0*0 / fax J ((vj) K - Vj) 

de z£K ^ ez 
n „ 



node zGA j=l 

For each node z G K, there exist a number M z < 1 and a sequence of elements K z m C 
K,m — 1, ■ ■ ■ , M z , such that iY 2 l = K, K zi and K Z) i + i have a common edge/face e^j, z = 
z,Mz '■= e z is an edge/face of K z>Mz . Clearly, we have 



1, ■ ■ ■ , M z — 1, and e 



\{ v j) K ~ v j \ ^ 


/ 










< 


/ 






/e z,M z 





e z,M z 



A' 



+ ••• + t 



Vj + 



1 



hz,l Je z ,i h Zj i Jr 

1 



- (*>i)jr., ) 



- (^')a z . m . ) 



< 



M z 

£ 



m=l J 9K Z , 



(«i> 



K ~ V 3 



- v, 



m=l 



L 2 (dK z , m ) 
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Thus, by the trace inequality and Poincare inequality, we obtain 

n M z 

-1/2 \^ / 1,-1/2 



IV < eh K 



node z£K j=l m=l 



+4 /2 



L 2 (K z , m ) 



3=1 



which, combining with (4.24)-(4.26), yields the result immediately. 



□ 



From Lemma 4J3 and Theorem 34^ we have the following result which gives H 1 ap- 
proximation estimates of the space Wh- The proof is omitted. 



Lemma 4.7. Let Uh := Z^ui. Then 
(4.27) 

(4.28) 



h 



- < I i l in I 1 / 2 1 I 

/i 3/2 , 



|«i - %lffi(n r , ) ~ e l u o|//2 (ni) + 



Wo 



where ttr h is defined in (2.19). 



5 Error estimates for the FE-MsFEM 

In this section we derive the if 1 -error estimate for the FE-MsFEM in the case where 
/3 = 1. For other cases such that /3 — 0, —1, the analysis is similar and is omitted here. 
Since the convergence analysis is only done for the periodic coefficient case, we will fix 
p = e in the later analysis. 

The following Lemma gives an inverse estimate for the function in space OMS(.K'). 

Lemma 5.1. Assume that vh G OMS(K). Then, we have 

- 1 + H- 1 ) 



(5.1) 



\ v h\h\k) < ( e 



Proof. Assume that vh = c&i ■ By the definition of ip" , we can extend the Vjj to the 
macro element S D K as following: vh = CiC^ipj, and HrVh = c i c ij<Pj, where c^,ipj, 



and ifj are defined in Section^ see (2.8 )— (2.10 ) j. It is easy to verify that vh satisfies 

-V ■ (a € (x)Vv H (x)) = 0, x E S, 

vh = T\. k vh-, x ^ dS. 



Hence, from Lemma 3.2 and dist(i^, dS) > Hk, it follows that 

\vh\hi(K) < (e _1 + H- 1 ) \\Vv H \\ L2{s) , 
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which yields 

\vh\h\k) < (e _1 + H- 1 ) \\VIi K v H \\ L2(s) < (e- 1 + H' 1 ) \\VYi K v H \\ L2{K) . 



Therefore, from Lemma 4.1, it follows the result (5.1) immediately. 



□ 



The following lemma gives the continuity and coercivity of the bilinear form A@(-,-) 
for the FE-MsFEM. 

Lemma 5.2. We have 



(5.2) 



\Af}(v,w)\ < 2 ||M|| 1Mr HMIIw* \/v,w e V KH - 



For any < 71 < 1, there exists a constant ciq independent of h, H , e, and the penalty 
parameters such that, if 70 > ato/ji, then 



(5.3) 



Ap(Vh,H,Vh,H) > \ \\\Vh,H\\\\hH V v h,H £ V h)H - 



Proof. (5.2) is a direct consequence of the definitions (2.3)-(2.6), (2.20), (2.21), and the 
Cauchy-Schwarz inequality. 

It remains to prove (|5.3l). We have, 



Ap{v h;H ,v h)H ) = Yl \\{^) 1/2 ^h,H 
KeM h , H 



.2 

\L^K)' 



v h,H 



2 J2 ■ n} [v hjH ] 

+ J2 (7 IIK^]lli2 (e) +7ie||[a e V^-n]||J 2(e) ) 

2 J2 Ji^Vv^H-n} [v h>H ] 



,h,H ^ 



— \\{a. e Vv hjH ■ n}\\ 2 L2 , e) 
7o 



It is obvious that, 



2 I {* € Vv h , H ■ n} [v KH ] < 2 \\{* e Vv h , H ■ n}|| L2(e) \\[v h>H )\\ 

edV h Je edT h 



L2( e ) 



eer h 



eer h 



Therefore, 



A/3(v h}H ,v htH ) > \\\v htH \\\l AH |e 



eer h 



L2( e ) 



It is clear that, for any e 6 T^, 
{a^v^a ■ n} 



(a e V^) ■ n + - [a e V^, H ■ n] 
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We have 

2 ^ — "TO 2 

(5.4) Ap(v htH ,v hiH ) > \\\v h ,H\\\ 1AH ~ l^y, IIK-ff]|| L2(e) 



eer h 



- ? E ( n( a£V ^) • + i n[ aev ^ • nlllL> 



L»(e) 



By the trace inequality, the inverse estimate (5.1), and e < H, we have 
(5.5) ]T || (a £ V^) • nf L2{e) = \\ (^Vv H ) ■ n\\ 2 L2(E) 



sS-E 



Cl ||/_ e Nl/2y- - 



e 



W(K E ) ' 



where Kg 6 is the element containing E. Therefore, from (5.4) and (5.5), it follows 
that 

2 ^ — To 2 

Ap(v h>H ,v h , H ) > \\\v h ,H\\\i A H ~ 2^ Ye H^]^ 2 (e) 

eer h 

- — ] M\l H - E 7ie ll[ afV ^ • n ]lli2 (e) 



7o 



27071 



> |||^,h||Ii feH - max 



eer h 
1 6C171 3 



2 ' 7o7i ' 27071 



\Vh,H\ 



i2 

\l,KH ■ 



Noting that 71 < 1, there exists a constant «o > independent of /i, if, e such that if 
7071 > ao then max ^ ? 27071 ) — I' This completes the proof of the Lemma. □ 

The following lemma is an analogue of the Strang's lemma for nonconforming finite 
element methods. 

Lemma 5.3. There exists a constant «o independent of e, h, H, and the penalty param- 
eters such that for < 71 < 1, 70 > aio/71, ^ e following error estimate holds: 



(5.6) 



K - U h,H\\\lh H 



< inf \\\u e - v h)H \\\ lhH + sup 



Vh,H£- V h,H 



Wh,H&V hyH 



\f n fwh,H dx - Ap(u e ,W h ,H)\ 
IH u, Mrllli,h,H 



Proof. For any v^jj G V/^h, from Lemma 5.2, (5.3), (2.14), and (5.2), we have, for < 
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71 < 1 and 70 > a /ji, 

\\\Uh,H - v h,H\\\l,h,H < Ap(Uh,H - V hiH ,U h ,H ~ V h)H ) 

= (f, u h ,H - Vh,n) - Ap(v h>H , u htH - v hjH ) 

= (/, u h ,H - v kh) - Ap(u e , u KH - v h>H ) + Ap{u e - v hjH , u KH - v hiH ) 

< Ui u h,H - V h , H ) - Afs(u € ,U h:H - V h>H ) + ||K ~ V h,H\Wl,h,H \\\ U h,H ~ V h>H \ \\ lhH 

Hence 

UK ~ U h,H\\\l,h,H - \W U £ ~ V h,H\Wl,h,H + \W U h,H ~ V h,Il\\\l,h,H 

I (/, U h ,H - V h>H ) - Ap(u e , U h)H - V KH ) I 



< \\W-v KH \h h „ + 



\U h ,H ~ Vh,H\\\lhH 



which yields the error estimate (5.6). This completes the proof. □ 

Now, we are ready to present the main result of the paper which gives the error estimate 
in the norm for the FE-MsFEM. 

Theorem 5.4. Assume that the penalty parameter < 71 < 1 and 70 > 0:0/71- Then the 
following error estimate holds: 



where fir is defined in (2.18). 



Remark 5.1. (a) The error bound consists of three parts: the first part of order 0[y/e 
+ H) from the oversampling MsFE approximation in f2 2 , the second part of order 

O ( — \ fl\ 1 1 ^ 2 ^) from the FE approximation in Qi, and the third part —=- — ^ + 

g2 ll^/llL 2 (Q r ) f rom the penalizations on T. 

(b) Suppose that the interface T is chosen such that dist(T, dfl) = 0(H). If the average 
value 

then we have 

^ r e u hHl/2 H 2 

\\K - Uh,H\\\l,h,H ^Vt+Jj+H+ —— + — . 

In this case, we may choose H ~ y/e and h ~ e 5 / 4 to ensure that \\\u e — u^hIWi k h ~ v^- 
The condition (5.7) may be checked by using the standard singularity decomposition results 
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for elliptic problems on polygonal domains IJlj [77| /. For example, we may show for the 
two dimensional case (n=2) that, if the inner angles of the polygon Q are less than |7r ; 
then (5.7) holds. 

Proof. According to Lemma |5.3 the proof is divided into two parts. The first part is 
devoted to estimating the interpolation error and the second part to estimating the non- 
conforming error. 

Part 1. Interpolation error estimate. We set vh,h as v^h^ = Uh,Vh,n\n 2 — ipH, where 
iih '■= Z^ui and ipH are defined in Lemma 4/T and Lemma 4^3 respectively. We are going 
to estimate |||it e — h H , i.e., to estimate each term in its definition (cf. (2.21)). 

First, from Lemmas 4.4, 4.7, we have 

1/2 



( 5 - 8 ) (iK - v h ,H\\l tH + ||(a e ) 1/2 V(M £ - Vh,H)\\ 2 L 2 (Ql) 

~ H |^o| ff 2(Q) + |wo|tyi,<x!(n) + jj l M olvyi,°o(n 2 ) + - |^i| 1 l' u olvKi,°°(n 1 ) • 
Further, since [u e ] = and [ux] = 0, it is easy to see that 



eer h 



eer h 



[Ui - u h 



(K) 



EaT H JE e6 r h Je 

By the trace inequality, we have 

/ (Ux ~4>h) 2 < H' 1 \\U X - tpH\\ 2 L 2 (K) + \\Ui-iPh\\l2(K) II v K -^h)\\ L 2 
Je 

Hence, taking a summation over Th, yields 

Y 7 / (Ml_ ^ )2 ~ He H Wl -^h\\ 2 L 2 {k) 

I / 2 \ 1/2 / 2 \ 1/2 

+ ~( Yl W U ± - ^v\\ 2 L2(K) ) ( W^( U ^-^h)\\ 2 L 2 (k) J ■ 



Thus, from from (3.7), (4.13), and (4.14), it follows that 

cv-r.. 6 JE 



JJ3 

— l«olfl3(n rH ) 
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Further, by the trace inequality, from Lemma 4.5 , it follows 



eer h 



J0_ 

e 



< 



< 



\ Yl ( k 1 IK ~ ^lli 2 (X e ) + IK - M\ L *{K e ) ll V ( M l - Uh)\\ L 2 iKe) ^ 



e&T h 



h 3 



\ u UH2(n rh ) > 



where flr h is defined in (2.19), and, in order to ensure the last inequality, for any vertex 
of elements in K^ h , we have chosen the corresponding edge/face in the definition of Scott- 
Zhang interpolation to be an edge/face of some element in K-p h - Therefore, from the 



above two estimates and Lemma 3.5, we have 

r3 



(5.9) 



— / [u e -V h>H \ < \Uo\ H 2 {Qr) + h \U \ 



2 



L 2 



Next, we estimate the term 



eer h 



7q 



[a £ V(it e - u/ifr) ■ n] 



(Or) 



\L 2 (e) 



It is easy to see that 



Y] — ||{a £ V(M e - v h H ) ■ n] 



By the trace inequality, we have 



lL 2 (e) 



7o 



— ||a e V(> e - u h ) • n| 



2 

L 2 (e) • 



|a e V(u e - Vh) • n||Ja (JS) < || V(u £ - «i) ■ n||^ 2(E) + || V(ui - ^) • n||J 2(E) 
< if -1 ||V(u e ~u x )f L2{K) + ||V(u 6 - «i)|| i2(J f) ||V 2 (m £ - ui) I 



+ /f- 1 ||v(«i - M\\ 2 L 2 (K) + ||v(«i - ^)|| L2( ^ ||v 2 ( Ml - 



(*0 
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Hence, a summation over Th follows that 



E — \\si e V{u £ - 4j h ) ■ n|| 2 5 



EdT H 



To 



<eH- 1 \\V(u e -u 1 )\\l H n TM ) + ^ 1 E ~ M 



\l 2 (K) 



KeK r , 



+ e ||V(u e - Mi)|| L 2 (QrH) ||V 2 (u e - ui) 

1/2 / 



^ 2 (f2r H ) 



+<{ E \Wui-m\\ 2 l2{k) ) ( e l|vv-v*C (J0 ) 

V K&K V ' ^ ifc tr„ / 



1/2 



Therefore, it follows that from Theorems |3.1[ [374] and Lemma [473 

e 



(5.10) 



- ||a e VK - M ■ n|£ 2(e) <(# + e) 2 Itiol 



To 



2 L. |2 



+ e 4 ||V/|| 



Similarly, by the trace inequality, we have 



E — l|a e V(u e - u h ) ■ n| 



eer h 



To 



2 

L 2 (e) 



£ e E ll a M«e - «i) • n||^ 2(E) + e ||a e V(wi - u h ) ■ n|| 2 2( , , 

IIV 2 ' 



Eer„ 



< eH- 1 1| VK - uOU^n^) + e ||V(u e - "OU^^ \\V 2 (u e - 
+ th~ 1 ||VK -M fc )||J 2(nrh) +e||V(wi -Uh)\\ L 2 (nVh) || VV|| L2(f7rJ 
Thus, from Theorems |3. 1|3.4 and Lemmas |3.5[ 4.7 it follows that 



(5.11) 



E- ||a e V(M e - ■ n| 



eer h 



To 

< f 2 17/ J 2 . 



2 

L 2 (e) 



6 \ u o\ H ^n) + e|w | + e 4 ||V/|| 

It is obvious that a same argument as above can be used to get the same error bound 
for the term 

E Tie || [a e V(w e - v KH ) ■ n]|| 2 2(e) . 



Thus, it follows from (5.8 )— (5.11 ) that 
(5.12) 



inf |||« e - v h>H \\\ lhH < H \u \ H2(n) + (v^+ jz) \u \ 

R 3/2 



h 



+ ^|fii| 1/2 kl 



+ 



—ff\ u o\m(n r )+ e l|V/|| ia( n r 
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Part 2. The non- conforming error estimate. Define 

Sjj := set of all interior edges/faces of M.h ■ 
For any Wh t H £ Vh,H, noticing that dVL 2 /T is empty, it is easy to see 

if,w h>H ) = (f,w h ) Ql + (f,w H ) n2 = (-V • (a e Vu e ),tw fc ) ni + (-V • (a e Vw e ), w H )n 2 
= (a £ VM £ , Vw h ) ni - J a e Vw e • nio,, 

+ V] ( (a £ Vtt £ , Vw H ) K - I a £ Vw e • n K w H J 
= (a £ Vu £ , Vw h ) Ul a ' VM ^ ' n ( w h ~ w H ) 

_r, J e 



Here the unit normal vector is oriented from K to K' and the jump [v] of v on an 
interior side E = dK D G is defined as [v] := v\k — v\k>- Furthermore, by the 
definition of Ap, it follows that 

Ap(u e ,w htH ) = {a e Vu e ,Vw h ) Ul a ' V ^ ' n ( w h ~ w h) + ^2 ( a ' V ^' ^ w h)k- 

eer h Je K£M H 

Thus, we have 

(/, Wh,H} ~ Ap(u e , w h , H ) = ~ ^2 a 'Vn e • n E [w H ] 

= - Y2 / a ' VM e • n E [w H - HhWh} 

= -^2 / af^Ue ■ n K (w H - R H w H ) - ^2 / ^u e -n(w H -U H w H ) 



Since 



k&Mh ujx egf h 
- Ri + R2. 

Ri = / f{w B - UhW H )- ^2 / ^^(^h - Rhw h ), 



we can estimate Ri by following j26j Chapter 6], or [151 Chapter 9], or the proof presented 
in [281 Theorem 3.1], and obtain 

(5.13) |Ri| < c(e\u \ H 2 {n , 2) + (Vi+ jj) |«olwn,« ( n a ) + e II/IIl^) ) 
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Next, we consider the second term R 2 . For any wh G Xh, it is easy to check that (see 
[151 [28]) 

w H = TIkWh + ex — t^t H m S, 

where Of G H 1 ^) is the boundary corrector given by 

-V.(a'Wf) = inS, «f | 8S = -x''^g^. 



By the Maximum Principle, we have 



l|0?IU»(s) < c||vn^H|| L ^ (5 ). 

Thus, from Lemma [4. 1[ we obtain 

\\w H - H-kWhWl^ok) < e||VIlxWH|U°°(5) < e|| VIIxWhIIl-^) 

< ei/- n / 2 ||Vn^|| L2(x) < eH- n ' 2 \\Vw H \\ L . {K) 

which yields 

l R 2 1 < Y \\ Vu ^\lhe) \\wh -H-hw h \\ L 2 {e) 
Ear H 

1/2 



£ l|VM e || L2(r) ( Rn W W H - ^HW H \\ LXJ{E) J 

( \ 1/2 



< eH 1/2 ||Vu e || L a (r) \\w H \\ l)H - 



On the other hand, by use of the trace inequality, it follows from Theorems 3.1, 3.4 and 
Lemma 13.31 that 

% Y W V ( U *- U ~I-)\\l2(K) + \\^( u e-Ul)\\ L 2 (K) \u e -Ui\ H2{K) 

k<lk Vh 

+ e 2 Y [ H X ! u o\h 2 (k) + \ u o\h 2 (k) \ u o\h3(k) J 

KeK rH 

2 2 3 1 1 . 2 

< e |m | H 2(q) + |uo|wi,cx>( n ) + e || V/|| i2 ( nr) • 

Therefore 

(5.14) |R 2 | < (e\u \ H2(n) + y/e\u \ whao(n) + e 2 ||V/|| L2(nr) ) H^H^. 
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It follows from (5.13) and (5.14) that the non-conforming error in Lemma 5.3 



sup 



fnf w h,H ~ Ap(u,w htH )\ 



u>h,H£Vh,H II \ w h,H\ || i t h,H 

£ e \ u o\ H i(n) + (y~ e + Jj) \ u o\wi,°°(n) + e 11/11^(0) + ^ II v /ll^(o r ) , 



which, combining with (5.12), (5.6), and (3.5), completes the proof. □ 



6 Numerical tests 

In this section, we first demonstrate the performance of the proposed FE-MsFEM by 



solving the model problem (1.1) with periodic and randomly generated coefficients re- 
spectively, and then show the ability of the FE-MsFEM to solve two multiscale elliptic 
problems with high-contrast channels. In all computations we do not assume that the 
diffusion coefficient values are available outside of the research domain. In order to il- 
lustrate the performance of our method, we also implement two other kinds of methods. 
The first is the standard MsFEM. The second one is a mixed basis MsFEM which use 
the oversampling multiscale basis inside the domain but away from the boundary, while 
use the standard MsFEM basis near the boundary. By this way, the mixed basis MsFEM 
doesn't need to use the outside information. 

For the methods FE-MsFEM and mixed basis MsFEM, the triangulation may be done 
by the following three steps. 

• First, we triangulate the domain Q with a coarse mesh whose mesh size H is much 
bigger than e. 

• Secondly, we choose the union of coarse-grid elements adjacent to the boundary dQ 
(and the channels if exist) as fli and denote \ Oi by For example, in our tests, 
we choose two layers of coarse-grid elements (and the coarse-grid elements containing 
the channels if exist) to form the domain Q\. Hence the distance of T away from dfl 
is 2H. 

• Finally, in Q 2 , we use the oversampling MsFEM basis on coarse-grid elements. While, 
in Qi we use the traditional linear FEM basis on a fine mesh for the FE-MsFEM, or 
use the standard MsFEM basis on coarse-grid elements for the mixed basis MsFEM. 
In our tests, we fix the mesh size of the fine mesh h = 1/ 1024 which is small enough 
to resolve the smallest scale of oscillations. 

Please see Fig. [T] for a sample triangulation. 

Since there are no exact solutions to the problems considered here, we will solve them 
on a very fine mesh with mesh size hf = 1/2048 by use of the traditional linear finite 
element method, and consider their numerical solutions as the "exact" solutions which 
are denoted as u e . Denoting by Uh the numerical solutions computed by the methods 



2G 



considered in this section, we measure the relative errors in the L 2 , L°° and energy norms 
as following 

\\ u h ~ ' u e\i,2 — U e \\ Lac \\Uh — U e \\ E 



L e\\]j2 \\"'e\\ijOo W^eWE 

In all tests, for simplicity, the penalty parameters in our FE-MsFEM are chosen as 70 = 20 
and 71 = 0.1. The coefficient a e is chosen as the form a e = a e I where a e is a scalar function 
and / is the 2 by 2 identity matrix. 



\u e 



6.1 Application to elliptic problems with highly oscillating coefficients 



We first consider the model problem (1.1) in the squared domain Vt = (0,1) x (0,1). 
Assume that / = 1 and the coefficient a e (xi,X2) has the following periodic form 



(6.1) 



a e (xi,x 2 ) 



2 + 1.8sin(2vrxi/e) 2 + 1.8 sin(27rx 2 /e) 



2 + 1.8 cos(27rx 2 /e) 2 + 1.8 sin(27rs 1 /e) 



where we fix e = 1/100. In our FE-MsFEM, we consider two choices of the parameter p. 
The first choice is p = e as stated in our theoretical analysis, while the other one p — h, 
the size of the fine mesh. The second choice is useful when the scales are non-separable. 
We first choose H = 1/32 and report the relative errors in the L 2 , L°° and energy norms 
in Table [T] We can see that the FE-MsFEMs give the most accurate results among the 



Table 1: Relative errors in the L 2 , L°° and energy norms for the model problem with periodic coefficient 



given by (6.1). e = 1/100, H = 1/32, h = 1/1024, 7o = 20, 71 = 0.1. 



Relative Error 


L 


L°° 


Energy norm 


MsFEM 


0.7263e-01 


0.7157e-01 


0.2560c-00 


Mixed basis MsFEM 


0.3422e-01 


0.3637e-01 


0.1714c-00 


FE-MsFEM p = e 


0.1238e-01 


0.1334e-01 


0.5159e-01 


FE-MsFEM p = h 


0.1252e-01 


0.1344c-01 


0.4840e-01 



methods considered here. Especially, when we take p = h, the FE-MsFEM still works 
well. 

The following numerical experiment is to show the coarse mesh size H plays a role as 



that describing in the Theorem 5.4 We fix h = 1/1024 and e = 1/100. Three kinds of 
coarse mesh size are chosen. The first one, H = 1/64, is denoted as 64 x 16; the second 
one, H = 1/32, is denoted as 32 x 32; the last one, H = 1/16, is denoted as 16 x 64. 
The results are shown in Table [2j From the table, it is easy to see that as H goes larger, 
the relative error in energy norm goes lower first and goes higher later, which is coincided 
with the theoretical results in Theorem 15.41 

Next we simulate the model problem with a random coefficient which is generated by 
using the random log- normal permeability field a e (x) by using the moving ellipse aver- 
age technique [20] with the variance of the logarithm of the permeability cr 2 = 1.5, and 
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Table 2: Relative errors of the FE-MsFEM with H = 1/64, 1/32, and 1/16, respectively, for the model 
problem with periodic coefficient given by (6.1 ). p = e = 1/100, h = 1/1024, 70 = 20, 71 = 0.1. 



Relative Error 


L' 2 


L°° 


Energy norm 


64 x 16 


0.1100e-01 


0.1502e-01 


0.9342e-01 


32 x 32 


0.1186c-01 


0.1290c-01 


0.5159e-01 


16 x 64 


0.1240c-01 


0.1795e-01 


0.6593c-01 



the correlation lengths h = h = 0.01 (isotropic heterogeneity) in x\ and X2 directions, 
respectively. The ratio of maximum to minimum of one realization of the resulting perme- 
ability field in our numerical experiments is 1.6137e+05. One realization of the resulting 
permeability field in our numerical experiments is depicted in Fig. [2j We also compare 




Figure 2: The random log-normal permeability field a € (a;). The ratio of maximum to minimum is 
1.6137e+05. 

three kinds of methods including the standard MsFEM, the Mixed basis MsFEM and the 
FE-MsFEM. In this test, we set H = 1/32 and p = h since there is no explicit e in this 
example. The relative errors for the three methods are listed in Table [3] From the table, 
we can also see that the FE-MsFEM gives the most accurate results among the methods 
considered here. 
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Table 3: Relative errors in the L 2 , L°° and energy norms for the model problem with random coefficient 
as shown in Figure^ H = 1/32, p = h = 1/1024, 7o = 20, 71 = 0.1. 



Relative Error 


L 2 


L°° 


Energy norm 


MsPEM 


0.3690e-00 


0.3731e-00 


0.6014e-00 


Mixed basis MsFEM 


0.1119e-00 


0.1857e-00 


0.4770e-00 


FE-MsFEM 


0.2635e-01 


0.8351e-01 


0.2975e-00 



6.2 Application to multiscale problems with high-contrast channels 

In this subsection, we use the introduced FE-MsFEM to solve two elliptic multiscale 
problems which have high-contrast channels inside the domain. 

In the first example, the coefficient a e is characterized by a fine and long-ranged high- 
permeability channel, which is set by the following way. The example utilizes the periodic 



coefficient a e in (6.1) as the background, while changing the values on a narrow and long 
channel that defined from (0.08, 0.49) to (0.92, 0.51) with new value a e = 10 5 (See Fig. g. 
For this problem, the "exact " solution is difficult to be obtained due to the singularities 




Figure 3: Permeability field 

near the corners of the high contrast channel. Our direct numerical simulation shows 
that the gradient values are large near the corners of the channel. We set H = 1/32 and 
p = e = 1/100. The results are presented in Table [4] where the relative errors in L 2 ,L°° 
norms as well as energy norm are shown. We observe that the FE-MsFEM performs 
better than other methods. 

In the second example, we use the coefficient depicted in Fig [4] that corresponds to 
a coefficient with background one and high permeability channels and inclusions with 
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Table 4: Relative errors for the model problem with the permeability depicted in Fig(3] p = e = 1/100, 
H = 1/32, h = 1/1024, 7o = 20, 71 = 0.1. 



Relative Error 


L 2 


L°° 


Energy norm 


MsFEM 


0.1640e-00 


0.2187e-00 


0.3773e-00 


Mixed basis MsFEM 


0.5415e-01 


0.2552e-00 


0.2977c-00 


FE-MsFEM 


0.1127e-01 


0.2090e-01 


0.6843e-01 



permeability values equal to 10 5 and 8 x 10 4 respectively. The results are listed in Table [5 

1 SS SS SS ss ss ss 1 ss ss is ss ss si ss ss ss l 
0.9 ■ ]] ® ® ® ™ ® ]] J' " S " S ™ ™ ]] J; 
0.8 ■ si s^^s^^^s^^^^^^^! s: 

0.5 - - ss ss ss ss ss ss ss ss ss ss ss ss ss ss ss 

0.4 ■ ™ SS 8 S S! " " " " " " S! " " " 

0.3 si 

0.2 ■ ™ 

0.1 - " ss " " ss " " " " " " ss " " " " 

I , , " , , 

0.2 0.4 0.6 0.8 1 

Figure 4: Permeability field: a e — 10 5 in two channels consisting of dark small rectangles; a e = 8 x 10 4 
in small square inclusions; a £ = 1 otherwise. 

We observe that our FE-MsFEM gives much better results than the other two methods. 
7 Conclusions 

In this paper, we have developed a new numerical scheme for the elliptic multiscale prob- 
lems which joints the oversampling MsFEM and the standard FEM together by using the 
penalty techniques. The idea is first to separate the research domain into two parts fl\ 
and f^2 = \ ^i that Qi contains the boundary dfl where the oversampling MsFEM can 
not apply, and singular points (or regions) where the oversampling MsFEM is inefficient. 
Then we apply the standard FEM on a fine mesh of Q\ and the oversampling MsFEM 
on a coarse mesh of SI2- The two methods are jointed on the interface T = dQi D d^2 of 
the fine and coarse meshes by penalizing the jumps of the function values as well as the 
fluxes of discrete solutions. 

A rigorous and careful analysis has been given for the elliptic equation with periodic 
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Table 5: Relative errors for the model problem with the permeability depicted in Fig|4j H = 1/32, 
p = h = 1/1024, 70 = 20, 7i = 0.1. 



Relative Error 


L 2 


L°° 


Energy norm 


MsFEM 


0.3546e-00 


0.4007e-00 


0.5943e-00 


Mixed basis MsFEM 


0.2243e-00 


0.2596e-00 


0.4997e-00 


FE-MsFEM 


0.4274e-02 


0.1284e-01 


0.7566e-01 



diffusion coefficient to show that, under some mild assumptions, if V is so chosen that 
dist(r, dQ) > H, then the if terror of our new method is of order 

0(^+^ + ^1 + 9. 

which exactly consists of the oversampling MsFE approximation error in Q2, the FE 
approximation error in Q\, and the error contributed by the penalizations on T. Note 
that, for simplicity, we have only analyzed the linear version of FEM for the discretization 
on Qi. 

Numerical experiments are carried out for the elliptic equations with periodic oscillating 
or random coefficients, as well as, the multiscale problems with high contrast channels, to 
verify the theoretical findings and compare the performance of our FE-MsFEM with the 
standard MsFEM and Mixed basis MsFEM. It is shown that, the FE-MsFEM performs 
better than the other two methods in all cases and much better in some experiments. 

There are several ways to improve further the performance of our FE-MsFEM. First, 
the linear FEM on Qi can be apparently extended to higher order FEMs to reduce the 
error term related to Secondly, since Qi may contains singularities, another interest- 
ing project is to consider a combination of adaptive FEM on local refined meshes on fl\ 
and oversampling MsFEM on Q 2 - Thirdly, based on existence numerical results for over- 
sampling MsFEMs [12], we conjecture that the theoretical assumption of dist(r, dQ) > H 
may be weaken to dist(r,<9f2) > Ce (at least, in practice) for some constant C. These 
will be left as future studies. 
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